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ABSTRACT ' ; 

The mid-span section of a previously reported 
control led-diffusion compressor stator has been experi- 
mentally evaluated -in cascade. Measurements were taken 
over a range of incidence angles for blade chord 
Reynolds numbers from 470 000 to 690 000. Blade chord 
length was 12.7 cm, aspect ratio was 2.0, and solidity 
was 1.67.’ Measurements included conventional cascade 
performance parameters as well as blade surface pres- 
sures. Computations were made for the inviscid flow 
field, ^surface boundary layers, and loss for several of 
the blade inlet 'angle conditions, and compared against 
corresponding data. 

NOMENCLATURE ■ : - . • ‘ ■ 

AVDR Axial Velocity-Density Ratio (Eq. 1) 

Cp surface pressure coefficient (Eq. 8) 

Cr local-to-reference pressure (Eq. 4) 

c' blade chord 

i ' incidence angle (air angle - blade metal angle) 

k local-to-reference mass flux (Eq. 3) 

P total pressure 

p static pressure 

q dynamic pressure 

Re Reynolds number based on chord 

s blade spacing 

V Velocity 

x/c fraction of chord 

B air angle 


y stagger angle (Fig. 5) 

p density 

a blade solidity . 

blade camber angle 
w loss coefficient (Eq. 2) 

Subscripts 

1 inlet plane 

2 outlet plane 

REF reference conditions 

t stagnation conditions 


INTRODUCTION .. 

In a paper presented at a recent International Gas- 
Turbine Conference, the design of a controlled diffusion 
compressor stator using an automated design procedure 
based on numerical optimization was described (1^). The 
controlled diffusion stator was a redesign of the first 
stage stator of the NASA Two-Stage Fan (2). The orig- 
inal design was highly successful, demonstrating a first 
stage peak adiabatic efficiency of 87 percent, and a 
radial distribution of loss across the stator which was.., 
remarkably low. The purpose of the redesign was to 
develop and demonstrate the feasibility of an automated 
design process using numerical optimization methods. 
Although significant performance improvement from an 
already excellent stator design was not expected, at 
least equivalent performance was expected. In addition, 
the redesigned stator was intended to provide a research 
vehicle for the new controlled diffusion class of. 

D lading. ... 

As described in Ref. 1, the controlled diffusion 
stator was designed using a series of computational 
analysis methods coupled by a numerical optimization 
procedure. The blade shapes were specified by a 
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geometry code extracted from the NASA throughflow com- 
pressor design program (3). Blade section mean lines 
and thickness distributions were described by several 
polynomial curves. The potential flow. about; each two 
dimensional blade section was calculated by the TSONIC 
code developed by Katsanis (4). Surface boundary 
layers were calculated using the McNally BLAYER code 
(5), and the optimization procedure used to couple all- 
o7 the aerodynamic and geometric codes together was the 
COPES/CONMIN code (6,7). 

The mid-span sectTon of the blade has been built 
and tested in a cascade wind-tunnel at the Naval 
Postgraduate School. Although inlet Mach numbers (0.2) 
in this subsonic tunnel are less than the design level 
(0.68), experiments can be conducted at realistic 
Reynolds numbers (500 000 to 700 000) on relatively 
large scale blading ^2.7 cm chord), on which detailed 
surface pressures can be obtained. Furthermore, since 
the computational codes used to design the blading were 
quasi-three-dimensional and steady state, it is of 
interest to obtain experimental data in a similar flow 
environment, as well as in the real, unsteady, fully 
three-dimensional environment of the compressor stage. 
The data obtained in cascade provide an excellent 
opportunity to assess the accuracy of the flow codes 
not only at the design point, but also at off-design 
conditions. In this report, data obtained over the 
full range of incidence angle (including flow visual- 
ization data), are compared against the calculated 
results from several flow codes. 

EXPERIMENTAL PROGRAM 

Description of the Test Facility 

The arrangement of the cascade wind tunnel facility 
is Shown in Fig. 1. The tunnel is supplied by a 700 HP 
olower in the basement of the building through an 
arrangement of sound baffles and turning vanes. The 
air enters the tunnel through two, two-dimensional bell- 
mouth contractions in series which produce a 152.4 by 
25.4 cm (60 by 10 in.) section just above floor level. 
Tne 25.4 cm width is maintained between heavy side 
walls, between which an adjustable test section is 
formed as shown in Fig. 2. The blades of the test 
cascade are mounted in a rack which can be moved along 
tne (fixed) back wall. The front wall is readily 
removed by overhead crane to give access to test 
section and instrumentation (Fig. 3). The air angle 
into the test cascade is governed by the setting of the 
adjustable lower end walls and the 59 inlet guide 
(or turning) vanes. The adjustable upper (exit) end 
walls are individually adjusted during operation to 
angles required to produce uniform* static pressure at 
the wall in the blade-to-blade direction. Static pres- 
sure taps are provided at 5.08 cm intervals in the 
Dlade-to-blade direction at vertical (axial) distances 
of 41.28 cm upstream and 16.51 cm downstream of the 
center of the test blades. Spanwise and blade-to-blade 
proDe traverses are provided at 31.12 cm vertically 
upstream and 27.31 cm vertically downstream of the same 
reference. Plexiglass windows in the removable wall at 
the test blading facilitate the use of flow visualiza- 
tion techniques (Fig. 4). 

Unusual features of the facility are the large 
scale, the large number of test blades (20, at 7.62 cm 
spacing) and the technique for controlling inlet air 
angle. Considerable facility development and test 
experience (8-13) preceded the present test program. 
Upstream flow uniformity, periodicity and steadiness, 
and acceptable axial velocity-density values without 
resort to boundary layer suction, were achieved as a 
result of the initial design features and subsequent 
experimental programs. Turbulence intensity, measured 


across five turning vane pitches at blade mid-height, 
was 2 percent. 

Test Blade 

The blade shape is represented in cascade in Fig. 5. 
The mean line is represented by two segments, each des- 
crioed oy polynomial curves. A symmetrical thickness 
distribution is built on the meanline. This thickness 
is made up of a front segment and a rear segment, each 
represented by a separate polynomial. The polynomial 
coefficients were arrived at through an optimization 
process described in Ref. 1. A summary of the section 
geometrical properties is given in Table I. 

Instrumentation and Data Acquisition 

The instrumentation is summarized in Table II. The 
flow survey probes were calibrated using Zebner's sur- 
face approximation method (Ref. 15) with computer pro- 
grams given in Ref. 16. 

A Hewlett-Packard HP3052 Data Acquisition system 
incorporating an HP9845A desktop computer and NPS/TPL 
HG-78K Scani valve Controller were used to acquire pres- 
sure and temperature data. Flow survey probes and tun- 
nel reference pressures were connected to one 48 port 
Scani valve and blade surface pressures and reference 
readings to another. Measurement uncertainties are 
given in Table III. Tunnel wall statics were displayed 
and observed on a 50 tube, water manometer. 

The china clay technique was used as described in 
Ref. 13. Surface patterns were recorded in color using 
a movie camera, and transition and separation distances 
were scaled from projected frames. The observations v' 
made with china clay were first 'calibrated.', by appli- . 
cations of the technique to both a sharp and a rounded 
leading edge flat plate model set into the cascade in 
place of the test blading (Ref. 13). 

Turbulence intensity measurements were taken with a 
full micron tungsten wire. The wire was calibrated 
using King's Law over fourteen points from a tunnel 
dynamic pressure of 0 to 39.4 cm of water. The corre- 
lation coefficient was 0.995. 

Testing Procedure and Data Reduction 

Seven inlet air flow angles from 24 to 46"" were 
selected and set in turn. At each angle setting similar 
procedures were followed. With the tunnel set to the 
desired inlet dynamic pressure (either 33 or 20 cm of 
water) tail boards and IGVs were adjusted to give near- 
uniform wall static pressures both upstream and down- 
stream of the test blading. Blade-to-blade probe sur- 
veys were then conducted over four adjacent blade spaces 
about the center blade. Probe data and reference data 
were recorded at 0.64 cm intervals over the outer 
spaces, and 0.25 cm intervals over the central spaces. 
The lower probe trailed the upper by 3.84 cm throughout 
the traverse procedure. A spanwise probe survey was 
then conducted at 1.27 cm intervals at the upper station 
with the probe located at about 2.5 cm from the suction 
side of the center blade. The survey procedure was 
repeated with the lower probe at approximately 2.5 cm 
from the pressure side of the central blade. The blade 
surface pressures and reference data were then recorded 
with tunnel reference data. 

China clay observations were made for similar set- 
tings of tne cascade at a dynamic pressure of approxi- 
mately 9.3 cm of water. The somewhat reduced velocity 
was necessary since higher velocities give too rapid a 
drying rate, and difficulty in recording and interpre- 
ting tne drying patterns. Blade static pressure distri- 
butions were recorded when china clay observations were 
made. These data could be compared with data already 
recorded at two higher Reynolds Number. 
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In the reduction of the probe survey data, after 
reduction to velocity, pressure and flow angle at each 
point, the data were referenced to quantities derived 
from tne tunnel plenum to atmospheric pressure (driving 
potential). As shown in Ref. 9 and 10, this permits 
integration of the measurements over the blade-to-blade 
displacement to obtain the Axial Velocity-Density Ratio 
(AVDR) and loss coefficient such that small variations 
in supply conditions during the survey have little 
effect. The AVDR was calculated using the definition 


ADVR 


£ <> 2^2 «2 

£ PlVi cos 8^ 


ds^ 


ds. 


and the loss coefficient using the definition 
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which is, equivalently. 
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Probe data were reduced using Eq.(3) to (6), using 
overlapping quadratic interpolation to carry out the 
numerical integrations to obtain AVDR and loss 
coefficient. 

Surface pressures were reduced to pressure 
coefficients using the definition 


Cp 


P - Pi 


- 

2 ^ 


( 8 ) 


To verify accuracy, the reference upstream conditions 
were derived in three different ways, namely: as the 

mass average of the blade-to-blade probe surveys; from 
the ensemble average of stagnation pressure readings 
and derived static pressure values during blade-to-blade 
prooe surveys; and as derived from the fixed Prandtl 
probe measurements. . 

CALCULATION METHODS 

Inviscid Flow Solution 

The inviscid flow about the blade section in the 
two-dimensional, blade-to-blade plane has been calcu- 
lated by two methods: the method described. by Katsanis, 
TSONIC (4), and the panel method developed by McFarland 
(j7). The TSONIC program solves the stream function 
equation by finite difference techniques for the sub- 
sonic, compressible flow regime. It is necessary to 
specify as input the fluid properties, inlet total tem- 
perature and density, weight flow, blade geometry, inlet 
and outlet flow angle, finite difference mesh, and a 
meridional distribution of streamtube height. .In the 
work presented herein, a linear distribution of stream- 
tube height was utilized so as to match the measured 
axial velocity-density ratio. 

Because the nature of the equations dictates that' 
the solution be of the boundary value type,' the outlet 
flow angle must be specified on the downstream boundary. 
In these calculations, measured exit angle was used. 

Tne TSONIC code uses a uniform mesh and, therefore, 
mesh packing in regions where more definition is desired 
is not permitted. To oetter define the calculated flow 
properties in the leading edge region, McFarland's panel 
code (^) was used. The code is a surface singularity 
method which solves the inviscid, irrotational , compres- 
sible blade-to-blade flow equations on a surface of 
revolution. Streamsheet thickness can be incorporated 
as a function of meridional distance. The governing 
equations are linearized by approximating compressi- 
bility effects, and solved using an integral technique 
(panel method). Up to 98 panel elements are possible. 
Computational time for a typical case, is under 4 CPU 
seconds on an IBM 3033 computer. The method is limited 
to subsonic flow and is less accurate for low solidity 
blade rows. The same input data set used for TSONIC 
can be used for the panel code. Because of user freedom 
to distribute panel elements, excellent definition in 
the leading edge region can be obtained. 

Boundary Layer and Loss Calculations 

Blade surface boundary layers were calculated using 
the program developed by McNally (5). In addition to 
the surface velocities, required input includes upstream 
flow conditions, fluid properties, and blade surface 
geometry. Among the output provided by the program are 
the conventional integral thicknesses, form factors, 
wall friction coefficient, and momentum thickness 
Reynolds number. 

The program uses integral methods to solve the two- 
dimensional, compressible laminar and turbulent boun- 
dary layer equations in an arbitrary pressure gradient. 
Cohen and Resnotko's method (^) is used for the lami- 
nar boundary layer, transition is predicted by the 
Schlichting-Ulrich-Granvi 1 le method (1^), and Sasman 
and Cresci's metnod (^) is used, for the turbulent 
boundary layer. 

A boundary layer which is initially laminar may 
proceed through normal transition to a turbulent boun- 
dary layer, or it may undergo some form of laminar 
separation before becoming turbulent. To provide 
f lexibi lity for analyzing this behavior, several 
program options are available to the user. The 
calculations may proceed from a laminar boundary layer 
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througn transition to a turbulent boundary layer. 
However, if laminar separation is predicted before 
transition, the turbulent calculations may be started 
oy specifying a factor by which the last calculated 
value of momentum thickness is multiplied (this factor 
is commonly chosen to be 1.0 to satisfy conservation of 
momentum). This new momentum thickness and a value 
for form factor based on the last calculated momentum 
thickness Reynolds number are used as initial values 
for the turbulent calculations. In no case is an 
initial turbulent boundary layer allowed which has a 
momentum thickness Reynolds number less than 320 (21) . 
Laminar separation is predicted when skin friction 
becomes zero. Prediction of turbulent separation is 
Inexact. A separation criterion common to compressor 
blade analyses which use integral boundary layer 
methods is the incompressible form factor. Hi. Values 
of 1.8 to 2.6 nave been proposed and used in the past 
(22) . In this report, critical values will be inferred 
from comparisons with data. 

Total pressure loss coefficients were calculated 
using Stewart's method (^) . Required input includes 
upstream flow conditions, fluid properties, displace- 
ment and momentum thickness at the trailing edge, blade 
spacing, trailing edge thickness, and exit velocity 
and flow angle. There is no satisfactory method for 
calculating loss when the boundary layer is separated, 
and a combination of analytic art and empiricism is 
generally used. In this report, for off-design cases 
where turbulent separation occurs, the displacement 
and momentum thicknesses at the location of critical 
form factor are assumed to remain constant to the blade 
trailing edge. These values are then used in the loss 
calculations from that point. Expected accuracy dimin- 
ishes as incidence angle increases or decreases 
significantly from the design point. 

RESULTS AND DISCUSSION 

Experimental Observations 

Since the cascade facility is somewhat unconven- 
tional, the quality of the flow is discussed before 
results are presented. 

Flow quality . Data were obtained at inlet air 
angles ranging from 25 to 46**. Care was taken at each 
setting to verify that the inlet pressure and velocity, 
and the outlet pressure were acceptably uniform, and 
that two-dimensional flow conditions with streamline 
\ contraction prevailed to the downstream measurement 
plane. Surveys of upstream and downstream flows for 
all test conditions are documented in detail in 
Ref. 14; In Figs. 6 and' 7 are shown, for illustration, 
the results of spanwise and blade-to-blade probe sur- 
veys, respectively, at the design air-inlet angle of 
39 . In Fig. 6 it can be seen that the flow remains 
nearly uniform over the center 40 percent of the blade. 

In Fig. 7, showing the blade-to-blade conditions, 
the variation which appears almost as noise on the 
inlet stagnation pressure profile is what remains of 
the wakes produced by the turning vanes. The peak-to- 
peak variation seen in Fig. 7 is about +1.4 percent of 
the reference dynamic pressure, or +1.1 percent of the 
cascade inlet dynamic pressure. Thus, the inlet velo- 
city was uniform to within 1.2 percent, and these small 
variations were accounted for by mass-averaging in 
evaluating both AVDR and losses. The inlet flow angle 
was uniform to within 0.5* over the central four blade 
passages which were surveyed. 

In viewing the downstream flow distributions in 
Fig. 7, it should be noted that data were taken at 
close intervals (0,25 cm) for integration over one 


blade passage, and over coarse intervals (0.62 cm) for 
the other three. All data points are shown connected 
oy straight lines. The minor departure from strictly 
repetitive conditions seen in the first (left-hand) 
blade interval was observed in all data sets. Inte- 
gration to obtain blading performance was therefore 
carried out over the third blade interval. Periodicity 
within the blade passages was verified using corre- 
sponding pressure taps on three instrumented blades. 

The quality of the test conditions did not change 
significantly except at the^two extreme angles tested. 
At the lowest air angle (25*), the inlet uniformity 
decreased fractionally as a result of operating the 
inlet turning vanes considerably off design, while the 
physical separation of the wakes was increased because 
of the less oblique wall angle. At the highest angle 
(46*), the uniform core of the spanwise profile was 
reduced from 40 to 20 percent of span as a result of 
increased wall and passage viscous effects associated 
with increased loading. 

Overall performance . The results obtained for the 
cascade performance based on survey probe measurements 
are shown in Fig. 8. Also shown for comparison are 
the corresponding curves for data obtained by Himes 
(12) for the equivalent cascade of DCA blades. Some 
aifferences are noted. First, the loss coefficient 
for the controlled diffusion blade was measured to be 
less than that for the DCA olade at the design point, 
and over the useful range of air inlet angle. The 
data for loss coefficient were well behaved and showed 
no discernible effect of Reynolds number, except 
(perhaps) at the highest air inlet angles. 

Tnere were measurable differences in behavior at 
air inlet angles less than design. The outlet air 
angle did not change significantly for the CD blade, 
and ranged from 1/2 to 3* greater than measured for 
the DCA blade. Larger values of AVDR were also meas- 
ured, although at angles approaching and exceeding the 
design value, the AVDR was measured to be lower than 
for the DCA cascade. 

The AVDR is largely a consequence of the behavior 
of the end wall boundary layers in passing through the 
cascade, and the reason for the high value at 33 . 
inlet air angle is not clear. However, a much reduced 
static pressure rise was also registered, and this was 
confirmedjby; manometer readings of the distributions 
of wall static pressures. In fact, the pressure rise 
across the cascade dropped very abruptly with reduced 
air inlet angle, the drop being associated with the 
aetection of leading edge separation on the pressure 
side of the blade. Subsequent examination of all the 
Information obtained in the tests -pointed to a need to 
obtain more information in the range of inlet flow 
angles between 30 and 36*. 

Vortex Shedding . When testing at the three lowest 
air inlet angles, pure, highly audible tones were heard 
as the tunnel was brought up to the required test con- 
ditions. The tones occurred at discrete settings of 
the plenum pressure, and the frequency increased (1080, 
1960, 2380, 2790, 3080 Hz) as the inlet velocity at 
which it occurred increased (approximately 30, 34.5, 

42, 48.5, 56 m/sec). Vortex shedding in the blunt 
blade wake was suspected when the Strouhal number based 
on trailing edge thickness was found to be 0.17 (when 
the tones occurred). An experiment was conducted to 
verify this explanation. The bluntness was removed by 
taping metal shim-stock to produce sharp trailing edges 
on all blades. The tones then did not occur. The 
removal of the shim-stock from one blade resulted in 
the reoccurrence of tones, but very faintly. 
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Vortex shedding in turbomachinery blade wakes has 
been observed or suspected in several turbomachinery 
research experiments to^). However, the occur- 
rence in the present case was particularly graphic 
Decause the frequencies were in the audible range and, 
despite the tunnel background noise level, some had the 
clarity of pure organ tones. 

In high speed compressors and turbines, these reso- 
nant frequencies can be over 100 kHz and therefore too 
high to be detected by other than the highest response 
transducers (see Ref. 28, for example). The issue of 
vortex shedding must be raised when attempting to obtain 
comparisons of data with analyses which aim at the cor- 
rect modeling of the flow. Whether any special treat- 
ment of the trailing edge flow is required when vortex 
shedding occurs, compared to when it does not, requires 
that a more detailed study be conducted than was pos- 
sible in the present series of tests. 

Comparison of Calculated and Experimental Results 

Inviscid flow computations, presented as static 
pressure coefficient versus percent of chord are com- 
pared with experiment in Fig. 9. The solid line repre- 
sents panel code calculations and the dashed line 
represents TSONIC results. 

Tnere is very little difference between TSONIC and 
panel code calculations, except in the trailing edge 
region. Even there, the differences are minor and are 
due to differences in the way each code models the 
trailing edge. No inviscid code will calculate the 
flow accurately near the trailing edge of a blade row 
because real, viscous effects are the most prominent in 
that region. For blades with round trailing edges, and 
particularly for the present case with a blade having a 
relatively large trailing edge diameter, the Kutta con- 
dition does not truly apply. Some artifice must be 
used, and code developers employ different methods. 

In TSONIC, Katsanis has used a mass injection or 
wake simulation model (1) sketched in Fig. 10. Tangents 
are formed at the intersection of the trailing edge 
circle and the blade surface, and extended to the ver- 
tical line which forms a tangent with the trailing edge 
circle (Fig. 10). The "wake" is then extended down- 
stream with an orientation determined by the downstream 
whirl boundary condition. Pressure is allowed to vary 
across the simulated wake. This modelling removes the 
trailing edge circle and the attendant accelerations 
and decelerations associated with curvature. The user 
may vary exit angle until "closure" of the pressure 
oistribution occurs at the trailing edge to simulate a 
Kutta condition, or tne user may specify an exit angle 
(such as an experimental ly measured angle) and accept 
the crossing of suction and pressure surface pressures 
near the trailing edge, as in the present calculations. 

McFarland sets the downstream static pressure by 
continuity considerations. He then forces the surface 
pressures on both surfaces to meet this pressure at the 
trailing edge by an interpolation process involving the 
last tnree panel elements on each surface. To achieve 
satisfactory results in the trailing edge region, it 
was found necessary in these calculations to use all 98 
panel elements, and to pack points in the trailing edge 
(FTE set equal to 1.05) . 

Both methods are artificial. But inviscid flow, in 
itself, is artificial, and comparisons will always dif- 
fer from the real flow in regions where viscous effects 
are strong. 

Experimental measurements are generally in good 
agreement witn both methods of calculation. Viscous 
effects near the trailing edge begin to appear (as a 
divergence between calculation and measurement) at the 
design point (Betal = 38.91) and increase as incidence 
increases. Tne last pressure tap on each surface is in 


a region of curvature and reflects the decreased pres- 
sure (increased velocity) there, particularly at lower 
inlet angles where suction surface boundary layers are 
relatively thin. Corresponding trailing edge accelera- 
tions do not appear in the calculations because the 
trailing edge has been modeled to remove them. 

Measured pressures compare very well with the cal- 
culations in the leading edge region. Because of mesh 
limitations, TSONIC did not pick up any leading edge 
acceleration/decelerations upstream of 3 percent chord, 
and therefore was quite limited in boundary layer cal- 
culations. The panel code provided quite adequate 
resolution. A point of concern to code users has always 
been the calculated leading edge pressure spikes, and 
they are evident here even at the design point (Betal = 
38.91). Tne excursions can be large and are mathematic- 
ally consistent. How much of such an excursion is 
actually felt by the real flow will determine the nature 
of the boundary layer behavior. It appears from the 
measured data that the flow does experience a large 
proportion of the calculated leading edge velocity dif- 
fusion. 

Two principal boundary layer calculations were made, 
one set using data from low Reynolds number tests to 
compare against flow visualization experiments performed 
at low Reynolds number (340 000), and one set using 
calculated pressure distributions corresponding to 
higher Reynolds number tests (670 000). The latter 
results were used to calculate loss coefficient to com- 
pare against measured loss. 

China clay flow visualization experiments were run 
on the suction surface of the blade over the full range 
of incidence angles, but at a low Reynolds number to . 
accommodate drying times of the solvent. The results 
are presented in Fig. 11. The axial position at which 
laminar separation-turbulent reattachment, transition, 
and turbulent separation is plotted against the inlet 
flow angle condition. The figure shows the existence 
of a laminar boundary layer to about 60 percent of chord 
at the two lowest flow angles. As incidence angle was 
increased, an abrupt change occurred and a laminar 
separation with turbulent reattachment was produced 
very near the leading edge (2 to 5 percent chord) at 
all of the higher inlet angles. 

The boundary layer calculations were made using 
experimental surface pressures as input, and results 
are also plotted on Fig. 11. In all cases, laminar 
separation with turbulent reattachment was predicted, 
rather than pure transition. The calculated laminar 
separation point is reasonably close to experiment for 
Betal = 28.00, but disagrees by about 15 percent chord 
at Betal = 32.87. The calculated trends appear to be 
more valid because as incidence angle is increased, the 
adverse pressure gradient on the suction surface becomes 
steeper. Since laminar separation is strongly influ- 
enced by pressure gradient, earlier separation would be 
expected at Betal = 32.87. At the higher inlet angles 
there is good agreement between experiment and calcula- 
tion. The adverse gradient near the leading edge is 
very steep and produces almost immediate laminar separ- 
ation. Due to the magnitude of the adverse pressure 
gradient in most cases, and the Reynolds number levels, 
the laminar separation bubbles are presumed to be short. 

The experimental static pressure distributions 
already reflect turbulent separation, i.e., reduced 
adverse pressure gradients in the separated region near 
the trailing edge. Using these reduced gradients as 
input to a boundary layer calculation is not likely to 
result in prediction of turbulent separation, and did 
not for these calculations. It should also be noted 
that tne experimental turbulent separation^ locations 
shown in Fig. 11 for Betal = 28, 32 and 36* are more 
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likely to be a localized drying phenomenon than a 
separation. 

Boundary layer calculations on both surfaces were 
conducted for the highest Reynolds number case over the 
full range of incidence angles using calculated surface 
pressures as input. The calculated laminar and turbulent 
separation locations are shown on Fig. 12, using the 
china clay results as reference. Predicted laminar 
separation locations occur earlier for both Betal = 28 
and 32.87, but agree with china clay results at the 
higner inlet angles. No turbulent separation 
is predicted for the two lowest inlet angles, but is 
predicted for all inlet angles 38** and greater. Assum- 
ing separation to occur at critical incompressible form 
factors of 2.0, 2.2, and 2.5, the corresponding loca- 
tions are plotted as percent chord in Fig. 12, The 
best correlation witn flow visualization data is given 
by Hi = 2.2. For the steep adverse gradient at Betal 
= 45.96, a flow visualization separation location of 
90 percent chord does not appear consistent. The 73 
percent chord value, corresponding to Hi = 2.2 agrees 
better with the flattening of the pressure distribution 
shown in Fig. 9(f). 

Total pressure loss coefficients were calculated 
using Stewart's method (^), and are compared with 
experimentally measured values in Fig. 13. Attempting 
to calculate loss at off-design conditions using an 
integral boundary layer method and a one-dimensional 
loss model is a dubious enterprise. The calculations 
obviously break down at the two extremes of inlet angle 
(solid line). But agreement is remarkably good at the 
other incidence angles. For cases in which turbulent 
separation is predicted, the values of displacement and 
momentum thickness corresponding to Hi = 2.2 were 
assumed to remain constant to the trailing edge of the 
blade where the loss model was applied. 

One notable point of disagreement is at the design 
inlet angle of 38.91. A recalculation of the boundary 
layer beginning at about 4 percent of chord beginning 
with a laminar boundary layer produces loss predictions 
in agreement with measurements (dashed line. Fig. D). 

To accomplish this it would be necessary for the boun- 
aary layer to have relaminarized in the reacceleration 
region between 4 and 25 percent of chord (Fig. 9(d)). 
Although no evidence of this is evident from the china 
clay experiments, this appears to be the only way the 
measured loss values could be reconciled with computa- 
tions. 

The boundary layer calculations for the lowest inlet 
angle (24.49) were particularly troublesome on the 
pressure surface due to the excessively steep adverse 
pressure gradient. Laminar separation followed by tur- 
bulent separation was predicted. The boundary layer 
calculations can be continued through modestly high 
form factors, and values of displacement and momentum 
thickness were obtained at the trailing edge. But, 
with such a steep pressure gradient, turbulent separa- 
tion would not be unexpected, and due to the much 
relaxed pressure gradient at 30 percent of chord, tur- 
bulent reattachment could be envisioned. Using this 
speculative approach, a loss coefficient comparable to 
experimental measurements was achieved by reattaching 
the turbulent boundary layer at 30 percent chord with a 
displacement and momentum thickness equivalent to a 
fourfold increase over values at predicted turbulent 
separation (dashed line. Fig. 13). Using such methods 
to preoict loss at far off-design conditions is clearly 
speculative. More sophisticated modeling is required 
or, at least, an extensive data base is necessary to 
permit consistent empirical enhancement. 


SUMMARY OF RESULTS 

The mid-span section of a controlled diffusion 
stator was tested in a two. dimensional cascade. 
Measurements over a wide range of incidence angles were 
obtained. A quite acceptable minimum loss coefficient 
level of 0.0241 was measured, a value lower than than 
measured on tne reference double circular arc blade. 
Operating range was slightly better than for the OCA 
blade. 

Measurements of surface pressure compare well 
against results computed by the TSONIC and panel codes 
over tne full incidence angle range. Integral boundary 
layer calculations using the McNally code showed gen- 
erally good agreement with flow visualization data when 
measured surface pressures were used as input. When 
computed inviscid surface pressures were used as input, 
laminar separation was predicted with less accuracy. 
Correlations between computed boundary layers and flow 
visualization data produced a critical incompressible 
form factor value for turbulent separation of 2.2 for 
this set of data. Loss calculations using the Stewart 
loss model provided remarkably good agreement except, 
as expected, at the extremes of the incidence angle 
range. 
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TABLE I - BLADE AND CASCADE GEOMETRY 


Camber, deg 41.09 

Max thickness, percent chord 7.0 

Leading edge radius, percent chord 0.9 

Trailing edge radius, percent chord 1.26 

Solidity . . . 1.67 

Stagger angle, deg 14.27 

Cascade aspect ratio 2.0 

Chord, cm 12.7 

Blade Spacing, cm 7.62 
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TABLE III - MEASUREMENT UNCERTAINTY 


Measurement 

Method 

Uncertainty 

Probe position 

Linear potentiometer 

±0.02 cm • 

Flow angle (Bi) 
(02) 

Angle potentiometer 
Angle potentiometer 

t.2 deg 
+.5 deg 

Pressures 
Plenum 
Wall static 
Blade surface 
Probe 

Scanivalve & transducer 
Scanivalve & transducer 
Scanivalve & transducer 
Scanivalve & transducer 

±.05 cm water 
±.05 cm water 
±.05 cm water 
±.05 cm water 

Flow field 
Stagnation pressure 
Static pressure 
Velocity 

Calibrated probe 

±.05 cm water 
±.2 cm water 
±.5 percent 





Figure 1. - Cascade wind tunnei test facility. 



Figure 2. -Test section instrumentation and physical dimensions. 








Figure 4 - View of test blading through plexiglas window. 
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gure 6. - Measured in 
in the spanwise dired 





BLADE-TO-BLADE DIRECTION, ciTi 

Figure 7. - Measured cascade inlet and outlet flow Held 
in blade-to-blade direction at design inlet angle. : 
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PANEL CODE (17) 

TSONIC(4) 

O EXPERIMENTAL, SUCTION SURFACE 

□ EXPERIMENTAL, PRESSURE SURFACE 



X/C, % chord 

(b) Inlet angle, 0^ = 28. 00 deg. 

Figure 9. - Comparison of experimental and 
calculated surface pressure coefficient. 
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PANEL CODE (17) 

TSON1CI4) 

O EXPERIMENTAL, SUCTION SURFACE 

-1.2 — □ EXPERIMENTAL, PRESSURE SURFACE 



X/C, % chord 

(d) Inlet angle, Pj = 38.91 deg. 
Figure 9. - Continued. 




(fhrilet angle, Pj = 45. 96 deg. 
Figure 9. - Concluded. 
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CHORD POSITION, X/C 


(a) Construction of effective wake; 



Figure 10. - Mass flow injection model: TSONIC program. 



AIR INLET ANGLf, deg 

Figure 11. - Comparison of boundary layer behavior between 
China clay experiments and calculations using experi- 
mental surface pressure distributions. Suction surface. 



AIR INLET ANGLE, pj, deg 


Figure 12. - Comparison of boundary layer behavior between 
China clay experiments and computations using calculated 
surface pressure distributions. Suction surface. 


CHINA CLAY _ 
(Re=3.4xl05) 
-- O — COMPUTED , 
(Re =7xl(F) 




Figure 13. - Calculated loss (Stew- 
art's model) cq^mpared wUh exper- 
imental measurements. 
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